Evidence for Unusual Dynamical Arrest Scenario in Short Ranged Colloidal Systems 
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Extensive molecular dynamics simulation studies of particles interacting via a short ranged at- 
tractive square-well (SW) potential are reported. The calculated loci of constant diffusion coefficient 
D in the temperature-packing fraction plane show a re-entrant behavior , i.e. an increase of diffusiv- 
ity on cooling, confirming an important part of the high volume-fraction dynamical-arrest scenario 
earlier predicted by theory for particles with short ranged potentials. The more efficient localiza- 
tion mechanism induced by the short range bonding provides, on average, additional free volume as 
compared to the hard-sphere case and results in faster dynamics. 
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Recently a large number of novel dynamical arrest phe- 
nomena have been described in systems where, besides 
the usual hard core, an attractive potential is present 
vifith a range much smaller than the hard core radius 
[Q-^ . This condition is not usually met in molecular liq- 
uid systems, where the attractive range is of the same 
order of magnitude as the hard core, but can be realized 
in colloidal systems where the size of the particles largely 
exceeds the range of the attractive part of the potential 
[Q-^ , and possibly many other systems of experimental 
interest, including globular proteins ||9|,[l0|]. From the ex- 
perimental point of view the most striking phenomenon 
associated with the attraction is the gradual disappear- 
ance of the liquid phase when the range of the attractive 
part of the potential diminishes. Ultimately, for narrower 
wells, the liquid-gas coexistence becomes metastable with 
respect to a crystal-fluid equilibrium, but nevertheless 
shows up as a metastable binodal curve in experiments 
and in simulations [Tl| ]. Various approximations to 
the liquid and crystalline free energies have been used in 
order to calculate the coexistence lines and for various 
forms of the attractive part of the interaction potential 

While the situation of the equilibrium phase diagram 
is quite clear, the metastable region of the phase-diagram 
related to the supercooled fluid is quite complex, and only 
partially clarified . It is now widely accepted that 

mode-coupling theory (MCT) provides a description of 
supercooled dynamics which is particularly suitable for 
understanding colloidal systems. MCT predicts in par- 
ticular the existence of a non-ergodicity transition which 
corresponds to a structural arrest, i.e. the impossibility 
of the particles of the system to move under the effect 
of the neighboring ones, the so-called cage effect p8| , p^ . 
Some measurements performed in colloidal systems sup- 
port the MCT predictions in more quantitative detail 



pO| , pl| . More recently MCT has been applied to systems 
interacting through a short ranged attractive potential, 
with the result of predicting a number of new and inter- 
esting phenomena pP, |2^ , p^ . Amongst the most strik- 
ing is the possibility to distinguish two types of transi- 
tion from a supercooled liquid to a glass, one mainly due 
to the repulsive part of the potential through the usual 
mechanism of the cage effect, the other to the attrac- 
tive part of the potential. The latter mechanism is due 
to the adhesiveness of the potential at short distances 
which produces 'clusters' of particles, a different mecha- 
nism of structural arrest than the 'blocking' or 'jamming' 
familiar in hard sphere systems. The attractive branch of 
glass transition curve that results extends from low val- 
ues of the volume fraction cj) to values of the order of the 
hard sphere transition and does not vary much with tem- 
perature. The repulsive glass curve passes from the hard 
sphere value at high temperatures to larger values of (f) as 
temperature decreases, thereby giving rise to a re-entrant 
phenomenon, i.e. the supercooled liquid phase extends 
into the glass region above the volume fractions of the 
pure hard spheres system. It is then possible, raising the 
temperature at constant volume fraction, to move from 
the attractive glass region to the supercooled fluid one 
and again in the amorphous hard-sphere-like glass region. 
We note in passing that the complementary possibility of 
driving the system across a glass transition both by in- 
creasing and decreasing the density has been shown to be 
possible in systems with long range interactions (Wigner 
glasses) 1 24 . In the case here under scrutiny, at high val- 



ues of (p, the two branches of the glass curves cross at an 
angle, and the attractive branch continues further into 
the glass region, giving rise to a remarkable coexistence 
of two different types of glass which can be characterized, 
for example, by their mechanical properties such as the 
shear modulus pH]. The glass-glass curve terminates in 
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a special singular point of the theory, named A3 in the 
MCT, where the relaxation processes have a peculiar be- 
havior |26j . Upon narrowing the width of the attractive 
potential, the glass-glass transition line tends to become 
shorter until it vanishes in a point which represents a high 
order singularity, an A4 singularity in the MCT language 
p2[ . The experimental and/or numerical verifications of 
these predictions are still scarce [§3[|j. A logarithmic 
decay, predicted by MCT close to the A3 point, was de- 
tected in particular in a micellar system at high packing 
fractions |2^. It is interesting to note the possible ap- 
plication of this type of results to the study of protein 
crystallization [islJgt]. 

We focus in this paper on a particular aspect of the 
attractive colloidal system, that is the extension of the 
glass transition line to high values of the packing fraction, 
in order to test one of the important predictions of the 
theory, i.e. the re-entrant behavior of the supercooled 
fluid-glass line. We simulate a monodisperse sample of 
N = 1237 particles of unit mass with a constant diam- 
eter (7 = 1 in a cubic box. The physical quantities are 
measured in units of the particle diameter a, the particle 
mass m, and the square well depth uq as unit of en- 
ergy. Temperature T is measured in units of energy, i.e. 
by setting the Boltzmann constant fcs = 1. With these 
choices, time is measured in units of a{m/uoY^'^ ■ The 
interparticle potential V{r) is the square well potential, 

V{r) — 00 r < a 

V{r) = -Uq a < r < a + a 

V{r) = r>a + A (1) 

The width of the attractive part of the potential is 
A = 0.0309, corresponding to a percentage variation 
e = A/(cr -I- A) = 0.03. We have investigated volume 
fractions from 0.10 to « 0.58 and temperatures from 
T = 0.32 to 50. For T lower than 0.32 the homoge- 
neous fluid phase is unstable with respect to gas-liquid 
phase-separation. MCT calculations, based both on the 
Percus-Yevick (PY) and Mean Spherical Approximation 
(MSA) structure factors ||2^, predict, for this specific po- 
tential, a reentrant fluid-glass line, as discussed in what 
follows. We have implemented the standard MD algo- 
rithm for particles interacting with SW potentials [p8| . 
Between collisions, particles move along straight lines 
with constant velocities. When the distance between the 
particles becomes equal to the distance for which V{r) 
has a discontinuity, the velocities of the interacting par- 
ticles instantaneously change. The algorithm calculates 
the shortest collision time in the system and propagate 
the trajectory from one collision to the next one. Calcula- 
tions of the next collision time are optimized by dividing 
the system in small subsystems, so that collision time 
are computed only between particles in the neighboring 
subsystems. 

For each run we have used particular care in equili- 



brating the system, and starting from an equilibrated 
configuration we have performed a simulation up to a 
time 10"^. Runs exhibiting the presence of crystalline nu- 
cleus at the end of the simulations were discarded. Since 
we simulate a monodisperse system, the formation of a 
crystal phase fixes the range of packing fractions where 
a stable (or metastable) fluid phase can be studied. 

We have calculated the self-diffusion coefficient D in 
the supercooled liquid phase via the long time limit of 
the mean squared displacement of the particles. For each 
of the 10 studied isotherms, D varies almost over three 
decades, showing a marked decrease at high volume frac- 
tions, before crystallization takes place. 

Fig.|l| shows the diffusion coefficient D as a. function 
of the packing fraction (f> = Nna^ / {6V) for the studied 
isotherms. The T = 50 isotherm diffusivity reproduces 
the hard sphere behavior. For each isotherm, simula- 
tions at larger volume fraction than the ones reported 
inevitably lead to crystallization during the run. Fig.^ 
shows that, when the kinetic energy is of the order of 
the potential depth (e.g. T = 0.75) , crystallization is 
shifted to larger packing fraction values as compared to 
both the hard-sphere case (e.g. T — 50) and to low T 
(e.g. T = 0.35). 

Fig.|^ shows the diffusivity behavior in the (f> region 
where the reentrant phenomenon takes place. Data are 
normalized by Do = ay^T/m, in order to take into ac- 
count the T dependence of the microscopic time. 

We see that on decreasing T at constant <j), D /Do first 
increases and then decreases again. Since the term 
in Do accounts already for the slowing down of the dy- 
namics associated to the different average particle veloc- 
ity, the increase oi D/Do on cooling must have a different 
origin. This peculiar feature can be explained as a re- 
sult of the competition between two different dynamical 
features produced by the increased bonding: (i) slowing 
down of the dynamics due to the formation of a larger 
number of bonded pairs and (ii) speeding up of the dy- 
namics due to the larger free volume resulting from the 
more efficient packing of the bonded particles — whose 
nearest neighbor distance is now imposed by the short 
range of the potential. The inset of Fig.^ shows all the 
studied state points. At all T, the low density limit co- 
incides with the hard sphere behavior. 

Fig. ^ reports the isodiffusivity curves, i.e. curves at 
constant D/Dq value, in the {<p,T) plane. It also re- 
ports the curve where crystallization takes place within 
the time of our calculation. The isodiffusivity curves 
shown in Fig. ^ can be considered as the precursor of 
the fluid-glass transition line that would take place when 
D/Dq ^ if crystallization would not occur first. The 
inset of Fig. || shows the MCT calculations for the same 
model as reported in [ p2[ . As in all cases studied previ- 
ously, the MCT calculation underestimates the location 
of the glass transition line but provides a correct frame 
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for understanding the reentrant behavior observed in the 
present calculations. 

In summary, we have established the presence of reen- 
trance in the dynamical arrest behavior of an extremely 
simple model, the square well one, when the width of the 
attractive potential is much smaller than the hard-core 
radius. This condition c;an be met in colloidal systems. 
This observed behavior, predicted by MCT, along with 
a number of other associated phenomena, is explained 
as competition between the hard-core caging char- 
acteristic of the hard-core systems — and the bonding 
caging which, in the case of very short range potential, 
localizes the particles in a more efficient manner than 
the hard core case. This stronger localization, imposed 
by incipient 'bonding' provides extra free volume which 
leads to more diffusional pathways. Further lowering of 
the temperature produces stronger cages, and the system 
then crosses over to the attractive glass scenario. It is an 
open challenge to find out in future if other MCT pre- 
dictions are supported by numerical investigation. This 
additional work, which requires the study of larger values 
of (j) and larger simulation time windows to better charac- 
terize the slow dynamics — based on SW binary mixture 
systems to prevent crystallization — is underway. 
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FIG. 1. Diffusion coefficient D as a function of the pack- 
ing fraction (f) for the e — 0.03 square well potential at several 
different T. 



FIG. 2. Scaled diffusivity D/Do (Do = a{T/m)°-'^) as a 
function of packing fraction </> for the SW potential for some of 
the studied isotherms. The inset shows an enlarged window 
to highlight the common low density limit. 
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FIG. 3. Iso diffusivity curves in the {4>,T) phase diagram 
with D/Do = 0.05(triangles), 0.04(squares), 0.025(crosses). 
The dashed line with filled circles represents the line where the 
system crystallizes within our maximum simulation time. The 
inset shows the theoretical MCT prediction for the ideal glass 
line (redrawn from Ref. [^) for both Percus Yevick (PY-fiUed 
diamonds) and Mean Spherical Approximation (MSA-open 
circles), separating the fluid phase from the glass phase. 
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